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ABSTRACT 

We perform a combined X-ray and strong lensing analysis of RX Jl 347. 5- 1145, one of the 
most luminous galaxy clusters at X-ray wavelengths. We show that evidence from strong 
lensing alone, based on published VLT and new HST data, strongly argues in favor of a 
complex structure. The analysis takes into account arc positions, shapes and orientations and 
is done thoroughly in the image plane. The cluster inner regions are well fitted by a bimodal 
mass distribution, with a total projected mass of Mtot = (9.9 ± 0.3) x 10^^ M^/h within a 
radius of 360 kpc/h (1.50. Such a complex structure could be a signature of a recent major 
merger as further supported by X-ray data. A temperature map of the cluster, based on deep 
Chandra observations, reveals a hot front located between the first main component and an 
X-ray emitting South Eastern sub-clump. The map also unveils a filament of cold gas in the 
innermost regions of the cluster, most probably a cooling wake caused by the motion of the 
cD inside the cool core region. A merger scenario in the plane of the sky between two dark 
matter sub-clumps is consistent with both our lensing and X-ray analyses, and can explain 
previous discrepancies with mass estimates based on the virial theorem. 

Key words: Galaxies: clusters: general - Galaxies: clusters: individual: RX J1347.5-1145 - 
X-rays: galaxies: clusters - Gravitational lensing 



1 INTRODUCTION 

Accurate determination of total mass of galaxy clusters is important 
to understand properties and evolution of these systems, as well as 
for many cosmological applications. Gravitational lensing, through 
multiple image systems (strong lensing) as well as from distortions 
of background sources (weak lensing), provides a reliable method 
to determine the cluster mass, which is inde pendent of the equilib- 
rium properties of the cluster (lMelliej[T99^ . The lensing mass de- 
termination can be compared to estimates based on measured X-ray 
surface brightness and temperature, w hich is instead based on the 
assumption of hydrostatic equilibrium (ISarazinlll988h or to dynam- 
ical estimates, which rely on the assumption of virialized systems. 
Combining optical. X-ray and radio observations of galaxy clusters 
is a major tool to investigate their intrinsic properties. In particular, 
the comparison of lensing and X-ray studies can give fundamental 
insigh ts on the dynamic a l state of the galaxy clu sters (see for ex- 
ample jAUenetaO l2002VD e Filippis et all (l2004l) ). on the vaUdity 
of the equilibriuni hypothes is and on their 3 -dimensiona l structure 
toe Filippis et al.Eo()5l : ISereno et al.ll2006l : ISerendl2007h . 



RX J1347.5-1145 (z = 0.451) is one of the most X-ray 
luminous and massive galaxy clust er known. This cluster has 
been the subject of numero us X-ray (Schindler etal.l ll99l 1 19971 : 
lAllen et al.l2002l : [Gittr&Schind ler 2004), op tical (Sahuetal. 1998) 
and Sunyaev-Zeldovi ch (SZE) effect studies (iKomatsu et al.ll200l] : 
iKitavama et al.ll2 004'). Formerly believed to be a well rel axed clus- 
ter, w i th a good ag reement between weak-lensing ( Fische r & TvsonI 
1 1991 iKling et al.ll2005r). strong lensing dSahu et al.lll998h and X- 
ray mass estimates dS chindler et al .1 [l997b , more recent investiga- 
tions revealed a more complex dynamical structure. In particular, 
a region of enhanced emission in the South-Eastern quadrant was 
first detected by SZE effect observations ( Komatsu et al. 2001) and 
later confirmed by X-ray observations that also me asured an hotter 
temperature for the excess component (lAllen et al.lE002) . This fea- 
ture has been inte r preted as an indication of a recent merger event 
(lAUen et al.l(200 3; uCitavama et al ]|20^. Furthermore, a spectro- 
scopic survey on the cluster members found a velocity dispersion 
of 910 ± 130 kms"\ which is significantly smaller than that de- 
rived from weak lensing, 1500 di 160 kms~^ (Fisher & Tyson 
1997), strong lensing, roughl y 1300 km and X-ray analyses 
1320 ± 100 km s"^ (aCohen & Kneibll200l see their table 4). A 
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Figure 1. VLT-FORSl image of the cluster RX J1347-1145. The black 
lines are the critical lines for the bimodal model (see Sect. 4) corresponding 
to the source redshift of Arc A4, z = 1.76. The white lines are the critical 
lines obtained for Arc A\ {z = 0.806). 



major merger in the plane of the sky was proposed as a lik ely sce- 
nario to reconcile all measurements ( C ohen & Kneibll2002b . 

In this paper, we further investigate the merger hypothesis 
by performing a combined strong lensing and X-ray analysis of 
archive data. We perform a strong lensing investigation based on a 
family of multiple arc candidates first proposed in Bradac et al. 
(2005) using deep VLT observations. Differently from previous 
studies, we take care of performing the statistical analysis in the 
lens plane, which is a more reliable approach than the source-plane 
investigation when working with only one multiple image system. 
We further refine our analysis by taking into account not only the 
image positions, but also the shape and orientation of the arcs. In 
addition, we exploit Chandra observations to gain additional in- 
sights into the dynamical status of the cluster, through spectral and 
morphological analyses of the X-ray halo, and to discriminate be- 
tween different evolutionary scenarios. 

The paper is organised as follows. Section [2 discusses the 
strong lensing image candidates selection from archive VLT and 
HST data. Section [3] describes the statistical method used in the 
lensing analysis whereas Sec. [His devoted to the X-ray data anal- 
ysis. Section [5] discusses the merger hypothesis. Summary and 
conclusions are presented in Sect. [6l Throughout this paper we 
use a flat model of universe with a cosmological constant with 

= 0.3 and Hq — 70 km s~^Mpc~^. This implies a linear 
scale of 5.77 kpc/^^ at the cluster redshift. 



2 OPTICAL DATA: CLUSTER MEMBERS AND ARC 
CANDIDATES 

The lensing analysis was based on archival VLT data, performed 
using the same dataset of Bradac et al. (2005), and more recent HST 
data. The VLT data were obtained with FORS 1 in high resolution 
mode (pixel scale 0.09'', total field of view 3.2' x 3.2'), using 
UBVRI Bessel filters. The seeing in the /-filte r is of 0.57" De- 
tails on the data reduction can be found in Brada5" etaD (l2005h . We 
generated object catalogs for each band, using SExtractor v2.4.4 to 
measure photometric and geometrical parameters. 

So far, more than five arc candidates for this cluster have been 
reported in the literature (see Fig[T])- The first two (Al and A4) 
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Figure 2. HST images of arcs A4, A5, AC (in the F814W band). On the 
left, we show the test curve A4; in the middle and right panel, we show 
the reference arc A5 and AC, respectively. The diamond corresponds to the 
observed arc A4 centroid. The crosses represent the counter-images of the 
four points sampling the test curve A4. 

were discovered by Schindler et al. (1995); later, HS T STIS im- 
ages revealed three additional ones, i.e. A2, A3, A5 JSahu et al.l 
1998). Recently, Bradac et al. (2005) reported the discover of sev- 
eral new arc candidates (AC, Bl B2, C, Dl, D2, D3, D4, E) using 
VLT data. Arc Al (2; = 0.806) is the only one with measured spec- 
troscopic redshift (Allen et al. 2002). Despite A3 could appear as a 
likely counter arc, this identification is still pretty uncertain. For the 
rest of the arcs only ph otometric redshift estimates are available so 
far (Br adac et aLl[2005 ). Thus, based on the existing literature we 
decided to start to use the most likely multiple-image system, e.g. 
the A4-A5 arcs, to model the cluster mass distribution. We then 
searched for additional images. Following Bradac et al. (2005), 
we compared photometric properties and used flux measurements 
in UBVRI to distinguish different arc families. We inspected the 
galaxies having I-band magnitude up to three magnitudes fainter 
than the BCG and with a distance to the BCG smaller than 1'. We 
confirm that the most plausible counter arc candidate is AC, giving 
a system of three images A4, A5 and AC (see Fig. (TJ which have 
consistent multi-band colors (using 5" aperture magnitudes) and 
surface brightness within the ei rors. Such a system is likely due to 
a source redshift at — 1.76 teradac etal ]|2003). It is important 
to stress that, while we can not exclude that there are additional 
multiple-image systems, our analysis indicates that there are no ad- 
ditional lensed images belonging to the A4-A5 system within the 
inspected region and above our detection limit. 

New observations from HST, in the F475, F814 and F850 
bands, became recently available (Proposal id. 10492). We anal- 
ysed the new data to check whether they confirmed out findings. 
The three images, obtained with the ACS camera, are not signifi- 
cantly deeper than the VLT data, but the higher spatial resolution 
(0.05" FWHM) allows to better disentangle nearby sources. The 
analysis confirms that A4, A5 and AC have consistent colors. In 
particular AC could now be resolved in a tangential arc and an 
overlapping source, see Fig. |2] We find that while the arc is fully 
compatible with the other images, the overlapping source is not, 
thus strengthening the results based on VLT data. 



3 STATISTICAL ANALYSIS 

We consider two density profiles to model the mass components: 
the first one for the cluster- sized halos and the second one for the 
galaxy- sized objects. We model the main mass (dark matter plus 
baryonic in tra-cluster ni edium) components as softened power law 
ellipsoids ( KeetonlEoOll) . The projected surface mass density at the 
projected position x in the plane of the sky is as usual expressed 
in terms of the convergence i.e the ratio of the local surface 
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mass density of the lens S(x) to the critical surface mass density 
Ecr = c^Ds/{4:7vGDdDds), being Ds, Dd and Dds the source, 
lens and the lens-source angular diameter distances, respectively. 
For a general softened power law ellipsoidal model, 

where b gives the mass normalisation, q is the projected axis ra- 
tio, a is the slope index, Oc is a core radius and ^ is an elliptical 
coordinate. We will adopt, for the main dark components, the non 
singular isothermal ellipsoid (NSIE) which has a fixed to 1 . 

We model the galaxy-size d (dark matter plus stellar content) 
halos as pseudo-Jaffe models (iKeetonI l200lh . The pseudo-Jaffe 
model is obtained by combining two softened isothermal ellipsoids 
in such a way that the projected density falls as outside the 
cut-off radius Ocut{> Oc). 

We performed a fit taking care of both the contribution 
from the image positions, Ximg^ and the arc shapes and orienta- 
tions, Xarcs^ SO th at = 'xf.rrr n + Xarcs- Wc cvaluatcd the m in the 
image plane ( Koc hane3ll99lh using the Gravlens cod^B ( KeetonI 
[2OOD. The ximg accounts for the agreement between the observed 
positions of the arc centroids and their predicted values. Comput- 
ing the Xhng in the image plane is the most reliable technique when 
dealing with just one multiple image system. In fact, the alternative 
procedure of estimating x^ in the source plane ( Kayser et al. 1990), 
despite of being computationally faster, does not take into account 
how many images make up the observed system, and might prefer 
in a flawed way mass models that yield a good fit adding fictitious 
images. 

We took into account the shape and orientation of the observe d 
arcs sampling each curve in a number of points (lKeeton|[200lh . 
A useful xiircs takes into account the distances between the fit- 
ted points of the counter-images of a test curve and the data points 
sampling the other reference curve. Such an approach is best suited 
to exploit the A4-A5-AC system, where the arcs are elongated and 
well defined. 

Our strong lensing statistical analysis was then implemented 
considering the system A4-A5-AC as a multiple image system. We 
considered both the centroid position and shape and orientation of 
the arcs as input data. We chose A4, the longest and better defined 
arc, as test curve, leaving the others as reference curves. The posi- 
tions of the sampling points are showed in Fig. [2 

The total number of constraints {N const) in our analysis is 
18 = 6 + 8 + 4, being 3 x 2 the centroid positions, 4 x 2 the 
additional sampled positions of the test arc and 2 x 2 the additional 
sampled positions of the Arc A5 . The image centroid of A4 and A5 
are used only for xlmg^ so that xlmg and x^arcs are independent. 
The number of parameters, Npar, and the number of degrees of 
freedom, Ndof — N const — Npar, for the different models are 
reported in Tab. [T] In our fitting analysis, we assume a positional 
uncertainty of 0.01'. 

3.1 One component 

As a first step, we considered the simplest tentative model, i.e a 
single main mass component centred in the neighbourhood of the 
BCG. Although we performed the whole investigation in the image 
plane, we first considered also the Xsrc minimisation in the source 

^ The software is available via the web 
site |http://cfa-www.harvard.ed'u7castles| 



plane. This allows i) to further check if other additional images of 
the chosen lens system are formed and ii) to allow an easy com- 
parison with previous works. 

A model with a single dark matter component provides a good 
fit only if two further images are present in the A4-A5-AC system. 
They appear when we check the Xsrc fitting result on the image 
plane. These two additional images should be lensed in the field of 
view of the VLT and HST observation, 0'.2 SW of the BCG, with 
amplification factors of the same order of the other three images 
and should then be easily detected, if present. This model has been 
then discarded. This first model is essentially the same as the one 
used in Bradac et al. (2005) as a first step in their fitting procedure. 

We then performed the x^ minimisation in the image plane let- 
ting all the dark component parameters vary ( Ndof — 10). A single 
mass component yields a very poor fit with a very high x^arcs- The 
mass component is located O'.IS East of the BCG. In Tab. [T]you 
can see the performance and the properties of this model, together 
with the other models discussed in the following. 

A single dark matter component model is therefore not appro- 
priate to describe the matter content in RX J1347.5-1145, hence 
revealing a possible more complex configuration. Even if the de- 
tails of the results in the following subsections strongly relies on 
the position of the counter-arc AC, the information contained in the 
shapes and orientations of A4 and A5 alone indicates a somewhat 
irregular mass distribution. 



3.2 Two components 

We tried to improve the fit by adding a second main mass com- 
ponent. We first considered the model with the two main matter 
components being fixed at the positions of the two brightest galax- 
ies ( Ndof — 8) - This model is similar to the one proposed by 
lAllen et al.l (l2002h . where the mass components were constrained 
by requiring that the overall potential was able to produc e the north- 
ern Arc 1. It is worth noticing that, as lAllen et al.l ( l2002h themselves 
remarked, their best model could not explain A4 and A5 as images 
of the same source, as instead suggested by later photometric ob- 
servations. 

When we consider the A4-A5-AC system as input data, the 
image positions are poorly reproduced and we get a radial A5 with 
an extremely high xiircs- Even if we move the position of the 
second mass component from the second brightest cluster galaxy 
(SCG) towards the peak of the south-eastern X-ray substructure, 
the fit is still very poor. 

We then made an additional step further, and let the position 
of the second dark matter halo free to vary, while leaving the main 
one fixed on the BCG ( Ndof = 6). A first acceptable result is that 
a second component at ^ 1^.2 East of the BCG can improve the fit 
substantially. The main improvement obtained with this model is 
that the Arc 5 is now tangential, even if too much elongated and 
shifted towards North with respect to what observed. Adding a sec- 
ond component clearly improves our fit, strongly suggesting a bi- 
modal structure for the cluster. 

We further explored the bimodal structure of the cluster by 
relaxing also the position of the first main mass component ( 
Ndof — 4), see Fig.|4] The fit obtained is pretty good, Xtot = 7.13. 
The two dark matter components have comparable masses and their 
centres do not coincide with either the BCG or SCG (see Fig.|4]). 
The first component is located North- West of the BCG, the second 
one South-East (but far away from the secondary X-ray peak). 



4 Miranda et al: RX J1347-1145 strong lensing and X-ray analysis 




Figure 3. The position of the source galaxy (triangles) of the arcs A4-A5- 
AC is shown together with the tangential caustic (dashed), radial (solid) 
curves and lips (dashed-dotted line). The right panel, where the naked cusp 
and the source are clearly visible, shows an enlargement of the left panel. 



3.3 Multiple components 

To test the reliability of the bimodal scenario and give an insight 
on some possible degeneracies, we explored some alternatives ac- 
counting for other physical effects. We then explored the effect 
of adding some further components, by considering a model with 
three main components, keeping the first two dark matter compo- 
nents fixed and adding a third free component, Ndof — 2. We 
fixed the first component in the BCG (coincident with the first X- 
ray peak) and the second one in the second X-ray peak leaving 
the third one free to vary. Even if the image positions are well re- 
produced, this model does not improve substantially the fit to the 
observations with respect to the model with two components. In 
fact, the mass clump centred on the X-ray peak is very small com- 
pared to the others and the global result is very similar to that with 
one component fixed and the second one free to vary. For all mod- 
els, when leaving only one component free to vary, while fixing the 
other/s one/s, the free component is always located in the eastern 
part of the system, showing a mass of the same order of the other 
main component, while the second frozen component, located ei- 
ther in the SCG or in the second X-ray peak, tends to be much 
smaller . 

We eventually considered the effect of galaxy- sized halos. 
Starting from our bimodal model with two free components, we 
added 10 galaxy-sized mass components centred at the location of 
the 10 brightest galaxies of the cluster (selected in the I-band). We 
expect that only those close to the images can affect the lens config- 
uration. We scale the p- Jaffe parameters w ith the luminosity using 
the following relations kneib et al.lll99i) : ao = cr*(L/L*)^/^ 
and rcut — r*^^^ (L/L*)^/^, respectively. We scale the core radius 
in the same way as Vcut. Note that all galaxies, except the BCG and 
SCG, are assumed to be spherical, see Tab. |2] The proportionality 
constants a* and Vcut are then fitted together with the parameters of 
the main mass components, so that Ndof — 2. The critical curves 
for this model are plotted in Fig. [T] (where we also show the critical 
curves corresponding to the Arc Al redshift, z — 0.806) and the 
caustics in the source plane are shown in Fig. [3] The best fit val- 
ues are listed in Tableland the surface mass density is plotted in 
Fig. 131 The errors on the parameters were estimated performing 100 
Monte Carlo simulations. Adding galaxy- sized halos causes a shift 
in the positions of the total mass distributions, towards the two main 
cluster galaxies. We note that when we introduce a galaxy-halo lo- 
cated at the SCG, the peak of the second main component moves 
South-East. Figure [2 reports three enlargements showing the data 
and fitted points for the system A4-A5-AC. 



3.4 Model comparison 

Lensing information strongly supports a bimodal structure. The re- 
liability of this result appeared by considering different alternative 
scenarios, which are summarised in Table [T] To further asses the 
goodness of our models, we exploited the Akaike (AIC) and the 
Bayesian information criteria (BIC) (see Liddle 2004 and refer- 
ences therein). AIC is defined as AIC = + '^Npar, while 
BIC — Npar log Nconst whcrc is the total for the 

model. The best model is the model which minimises the AIC or 
the BIC. By looking at Tab. \T\ we can see that the one component 
model is highly disfavoured, while adding a second mass compo- 
nent improves the fit, since the BIC value substantially drops down, 
suggesting a bimodal structure for the cluster. When we let one 
component position vary, the BIC value becomes even smaller and 
a substantial improvement is obtained when you let both compo- 
nent positions be free. Adding the 10 galaxies does not improve 
considerably the fit, but it helps to better refine the features of 
the cluster model. In other words, a bimodal model is naturally 
favoured by our analysis, which definitely rules out the case with 
just one component. 

The consideration of different lens configurations allowed also 
an insight into model degeneracies. As can be inferred by the pre- 
vious analysis, the main degeneracy is connected to the positions 
of the two cluster- sized halos making up the binary lens. The fit 
does not change in a dramatic way if the two centres are translated 
keeping the separation fixed at ~ l^ 

Our result on the model configuration does not depend on the 
estimated source redshift. In fact, a variation of Zs implies a rescal- 
ing of the critical density, but does not affect the positions and the 
relative mass ratios of the components of the cluster. We checked 
that our result does not change in a significant way if we employ a 
different sampling scheme for the curves. In the same way, adopt- 
ing either larger or smaller positional errors does not alter the fit 
values. 



4 X-RAY ANALYSIS 

We looked for further evidence of the complex structure of 
RX J1347.5-1145 by analysing the intra-cluster gas at X-ray wave- 
lengths. RX J1347.5-1145 was first detected in the X-rays by 
ROSAT in the All-Sky Survey, which showed and extremely X-ray 
bright extended source ("Schindler et al."l995'). Subsequent ROSAT 
HRI and ASCA pointed observations revealed it as the most lu- 
minous X-ray cluster known at that time, with a bolometr ic lumi- 
nosity Lx ^ 2 X 10^^ ergss"^ dSchindler et al.lll997h . Based 
on a ~ 20 ksec Chandra observations, Allen et al. (2002) re- 
ported the discovery of a region of relatively hot, enhanced X-ray 
emission, approximately 20^^ to the south-east of the main X-ray 
peak (which is located at RA= 13^47^30^6, Dec= -11°45'09'.'3 
(J2000)), separated from the main cluster core by a region of re- 
duced emission. Later XMM-Newton observations confirmed these 
results dGitti & Schindledi2004l : lGitti et aLll2007h . 

In this work we analysed an archival Chandra observation of 
RX J1347.5-1145, carried out on September, 2003, using ACIS-I 
(already published by Allen et al. 2004 within a study of relaxed 
galaxy clusters). The data reduction was performed using CIAO 
3.3.0.1. The net good exposure time, after removing all periods 
of high background, is of 56 ksec. This longer Chandra exposure 
roughly confirms previous findings detailed above. In the following 
we will hence exclusively point out some new results, compared to 
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Model 


Npar 


Ndof 




Xarc 


^^tot 


AIC 


BIC 


One Comp 


8 


10 


244 


2827 


3061.00 


3077.00 


3084.12 


Two Comp: 2 fixed 


10 


8 


228 


770 


998 


1018 


1026.93 


Two Comp:l free 1 fixed 


12 


6 


39.52 


89.15 


128.67 


152.67 


163.35 


Three Comp: 2 fixed 1 free 


16 


2 


35.36 


72.84 


108.20 


140.2 


154.44 


Two Comp: 2 free 


14 


4 


0.48 


6.65 


7.13 


35.13 


47.59 


Two Comp: 2 free +10 galaxies 


16 


2 


0.13 


6.03 


6.16 


38.16 


52.41 



Table 1. Summary of our lensing models. The values, the Akaike (AIC) and Bayesian (BIC) information criteria are reported. Note that Npar counts 2 



free parameters due to the unknown source position during the fitting procedure. 



Cluster-Size 


b (arcmin) 


xi (arcmin) 


X2 (arcmin) 


e 


PA (deg) 


(arcmin) 


Ocut 


IDMcomp 


0.872 ± 0.03 


0.333 ±0.04 


0.522 ±0.03 


0.08 ±0.04 


-12.4 ±0.2 


0.521 ±0.03 




IDMcomp 


0.747 ±0.03 


-0.494 ±0.02 


-0.277 ±0.01 


0.05 ±0.03 


27.34 ± 0.09 


0.40 ±0.03 




Galaxy-Size 


cr (km/s) 










^c(kpc) 


Ocut(kpc) 


BCG 


250 ± 100 


[0.0] 


[0.0] 


[0.3] 


[-0.10] 


[0.5] 


100 ± 20 


SCG 


212.5 ± 100 


[-0.301] 


[-0.033] 


[0.23] 


[33] 


[0.425] 


85 ±20 



Table 2. Parameters for the dark matter halos and of the two brightest galaxies of our best fit model. Values in square brackets are fixed and not minimised, e 
is the ellipticity; PA the position angle measured North over East. The coordinate system is centred on the BCG galaxy. 



what already published in literature. For the extraction of the az- 
imuthally averaged temperature profile (see ^4.1K obsID 506 and 
507, with net exposure times of 9 and 10 l:sec, respectively, were 
also used. In the innermost regions of the cluster the hot gas closely 
follows the BCG gravitational potential, with the X-ray emission 
approximately oriented in the NS direction (Fig. [5]). The main X- 
ray peak is slightly shifted towards SW respect to the centroid of 
the overall cluster emission. 



4.1 X-ray Spectroscopy 

Throughout this section, spectra have been modelled (in the 0.3 — 
7.0 keV band) to an absorbed XSPEC isothermal plasma emis- 
sion code bv Kaastra & Mewe (1993), including the FeL calcu- 
lations of iLiedahl et all (119951) (IMekaL) folded through the ap- 
propriate response matrices and corrected for the ACIS time de- 
pendent absorption due to molecular contamination. The absorb- 
ing column density was fixe d at the nominal Ga lactic value uh = 
4.85 X 10^° cm"2 (Dickev & Lockman|[T99Qb . All errors, if not 



otherwise stated, are 90% confidence limits. Within 1.5^ from the 
main X-ray peak (to the full 360°), we measure a mean emission- 
weighted temperature of kT — 11.7 ± 0.5 keV and a metal abun- 
dan ce Z — 0.30 ± 0.06 Zq in a^eement with prev ious Chan- 
dra dAUen et al.ll2002h and ASCA JSchindler et al.|[l997.) measure- 
ments. 

The azimuthally averaged temperature behaviour was in- 
spected by extracting spectra into elliptical annular bins centred 
at the cluster X-ray main peak, and excluding the South-Eastern 
quadrant. We observe a sharp temperature decrease within the cen- 
tral ~ 200 kpc, consistent with the presence of a massive cool 
core, and a more gentle decline in the outer regions d Allen et al.l 




100 1000 
r [kpc] 



Figure 6. Projected X-ray gas temperature profile in the 0.3 — 7.0 keV 
energy range. 
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Figure 4. Projected mass distribution as inferred from the strong lensing analysis of the system A4-A5-AC for some of the models discussed in the text. Left 
panel : multiple smooth dark components, one component fixed in the BCG, the second one fixed in the second X-ray peak, and the third one being free to 
vary. Middle panel: two smooth and free dark matter components. Right panel: our best fit model, composed by two smooth dark matter components plus the 
10 brightest galaxies. The coordinate system is centred at the BCG; units are in arcminutes. The surface mass density is in units of Ecr = 3100 /iM© pc~^; 
contours represent linearly spaced values of the convergence k from 1.3 to 0.5 with a step An =0.1. The two star shaped symbols denote the positions of the 
BCG and the SBG; filled dots represent the centres of cluster-sized mass components; the position of the SE X-ray sub-clump is given by an empty box. 




Figure 5. VLT-FORSl image of RX J 1347+ 1145. Superposed are the contours representing the lensing mass distribution (left panel) and the X-ray surface 
brightness (right panel). The units are the same as in Fig.l 



l200llGitti & Schindlei][2004 l) (see Fig.[6]). The deprojected tem- 
perature profile follows the same trend, and does not differ signifi- 
cantly from the projected one shown in the plot. 

We then performed a two-dimensional spectral analysis of 
the central 1.6^ x 1.6^ of the cluster. This was done with 
the aid of t he WVT binning algorithm by S. Dighl and T.S. 
Statler ( Diehl & Statlei]l2006l : Eappellari & Copinll2003h . and then 
extracting spectra for each resulting Voronoi cell (each containing 
at least 900 photons (S/N-^ 30) in the 0.3 - 7.0 keV energy range, 
after background subtraction). The resulting temperature map is 
shown in the left panel in Fig.[7l 

We also produced an additional temperature map using a to- 
tally different b inning technique: the contour binning algorithm by 
ISandersI 1200^ . This method chooses regions by following con- 
tours on a smoothed image of the cluster; the generated bins closely 
follow the surface brightness distribution. As above, spectra were 
extracted from each resulting region, each containing at least 900 



photons The resulting temperature map is shown in the right panel 
in Fig. [T] The relative errors of the temperatures in two resulting 
temperature maps are of the order of 10 — 20%, with a slight de- 
pendence on the temperature. Regions with larger errors were cut 
off from the final maps. 

The two temperature maps reveal a cool core {kT ~ 5.8 keV) 
which appears slightly flattened, extending from North to South 
(see left panel in Fig. [7]). Out to r ^ 1^ the temperature is approxi- 
mately constant around ~ 14.0 — 17.0 keV. Several much hotter re- 
gions are observed, one of which lies surprisingly close to the cool 
cluster core (the ''hotfronf at 9" south-east of the cluster centre). 
This hot front is located in the region of reduced X-ray emission 
that separates the cluster core from the the X-ray sub-clump, which 
instead shows a slightly lower temper ature. The X-ray s ub-clump, 
contrary to what previously thought (lAUen etaLll2002h . does not 
therefore correspond to the temperature peak in the region. In the 
hot front region, neither optical nor X-ray point sources are ob- 
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Figure 7. Adaptively binned temperature maps of the central (1.6^ x 1.6^) performed using a Voronoi binning algorithm (left panel) and a contour binning 
algorithm (right panel). The temperature scale on x-axis is expressed in keV. Superposed are the logarithmically spaced smoothed X-ray contours. 



served. The same features are observed in the temperature map 
produced with the contour binning technique (right panel in Fig.|7]), 
showing that the result cannot be ascribed to the binning technique. 
A trail of cold gas (^ 8.0 — 10.0 keV), approximately 130 kpc long 
and extending form the BCG toward N, is also visible. More than 
half of the detected photons at that position will be projected from 
the surrounding gas; the intrinsic temperature of the cold trail is 
hence likely to be lower than the values quoted above. The cold fil- 
ament is probably a cooling wake, caused by the motion of the cD 
galax ies inside the cool core region (jPavid et al.l 19941 : [pabian et aP 
1200 ih. 



4.2 X-ray and SZE face to face 

Our X-ray analysis has revealed hot gas {kT ^20 keV) in the SE 
quadrant, associated to a hot front located between the cool core 
and the SE sub-clump rather than to the sub-clump itself. This front 
is most probably what remains of a shock front, a clear indication 
of a recent merger event. 

SZE effect studies, despite their low spatial resolution, further 
strengthen the above conclusions. The SE sub-clump of en- 
hanced emission was first detected through SZE measurements 
by Komatsu et al. (2001) . Subsequent X-ray and SZE ob- 
serva tions revealed a sub stantial agreement with previous SZE 
data (iKitavama et al.ll2004 and references therein). 
In order to characterise the high temperature substructure, Ki- 
tayama et al. (2004) performed a combined X-ray-SZE analysis. 
They assumed that the hot substructure was embedded in an ambi- 
ent gas identical to that in the other directions, deriving a temper- 
ature of 28 =b 7 keV for the hot component, in agreement with 
our X-ray measurement for the hot front. 



5 A MERGER HYPOTHESIS 

Our mass estimates for the inner regions in RX J1347+1145 are 
in substantial agreement with previous lensing analyses. From our 
bootstrap error analysis, we estimated an integrated total projected 
mass of (1.42 zb 0.05) X lO^^Mo inside a radius of 1.5 arcmm. 



consistent with the estimate of (1.2 it 0.3) x 10^^ from a com- 
bined weak and strong lensing analyses in Bradac et al. (2005). We 
should point out that our results depend upon the correct redshift 
determination and identification of members of the multiple image 
system we use. An error of ±0.5 on the source redshift implies an 
error of ±1.0 x lO^^M© on the mass estimate. 

A first weak lensing result by Fisher & Tyson (1997) gave a 
total integrated mass of (1.5 ± 0.3) X lO^^Mo within - b' and 
a corresponding velocity dispersion of ctv = 1500 ± 160 km 
assuming singular isothermal sphere modefl This result was later 
confirmed considering a larger field of view by Kling et al. (2005), 
which obtained ctv = 14001}4q under the same assump- 

tions, while our mass estimate translates into a velocity dispersion 
cTv = 1620 ± 30 km s~^. Results from both our lensing and X-ray 
analysis indicate that the cluster is not in virial equilibrium; the in- 
ferred estimate of ctv has hence been derived only for comparison 
with weak lensing results. 

Since our analysis showed indications of dynamical complex- 
ity, estimates of the cluster mass as inferred from X-ray data under 
the hypothesis of hydrostatic equilibrium should be treated with 
caution. Nevertheless, just for comparison with the lensing mass 
we performed an X-ray analysis excising from the X-ray data the 
perturbed South-East quadrant , using both parametrical (single and 
double (3 models) and model independent approaches. We obtained 
an X-ray estimate of the mass within 1.5 arcmin of (8 ± 1 ± 2) x 
lO^^Mo, where the first error is statistical and the second one ac- 
counts for the uncertainty on the intrinsic geometry of the cluster, 
which might be either oblate or prolate. The total 3 -dimensional 
mass within a spherical region of radius 1 Mpc (2.9 arcmin) is 
Mtot = (1.1 ± 1 ± 2) X lO^^M©. The se estimates are in good 



agreement with previous XMM- Newton dOit ti & Schindlej 
Gitti et al.l l2007h and Chandra (lAUen et all 2002. 20041. 



2004 



2007 



Schmidt & Allenll2007h X-ray analyses, and are consistent as well 
with our results from gravitational lensing. However, X-ray mass 
estimates of perturbed system may be strongly biased. On one 
hand, neglecting the X-ray South-East sub-clump implies that we 

^ The values for fjv and M have been converted to our reference ACDM 
model. 
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are not considering an additional contribution to the total mass, 
which is consequently underestimated. On the other hand, the tem- 
perature in a merging system might be inflated with respect to the 
temperature of a relaxed cluster, which implies an over-estimate of 
the mass. 

Whereas X-ray and lensing estimates are in good agree- 
ment, the dynamical mass estimate is considerably smaller 
dCohen & Kneibll20 02Y Cohen & Kneib (2002) suggested a merger 
scenario to solve this puzzle. Until the merger is complete and 
galaxy orbits have virialized to the new total mass, the dynamical 
mass estimate would be biased toward the mass of the larger clump. 
Then, a major merger between two clumps of comparable mass 
could reconcile the discrepant observations. The bimodal mass dis- 
tribution inferred in our analysis is in agreement with this merging 
scenario . Moreover the merger is possibly happening in the plane 
of the sky, since the redshift distribution of the cluster members 
does not show any feature along the line of sight. 

The merger scenario is also supported by the X-ray and SZE 
observations of the intra-cluster medium. Both the gas sub-clump 
and the hot front in the SE quadrant are indications of a merger pro- 
cess. The overall agreement between the X-ray and lensing mass 
estimates might suggest that the emitting gas has had time to viri- 
alize. Hydrostatic equilibrium is also supported by the fact that 
the overall mass, temperature and luminosity of RX J 1347-hll45 
relate as expected by a typical for massive cluster (I Allen et al.l 
I2OO2'). A major merger scenario was also considered in Kitayama 
et al.(2004), which discussed as an head-on collision of two ^ 
5 X 10^^ M© clusters with relative velocity of ~ 4000 km 
would result in a bounce-shock with kT ^ 25 keV, as detected in 
the SE quadrant. 



6 SUMMARY AND CONCLUSIONS 

We have analysed both the lensing and X-ray properties of 
RX J1347. 5-11.45, one of the most luminous and massive X-ray 
cluster known. Based on the analysis of an arc family, photometri- 
cally selected, we have estimated the total cluster mass distribution, 
within a radius of ~ 500 kpc from the cluster centre. We performed 
a analysis in the lens plane and scrupulously modelled the arc 
configuration. A model with two smooth dark matter components 
of similar mass accurately reproduces the observations and yields 
a mass estimate in agreement with previous strong and weak lens- 
ing and X-ray studies. Our strong lensing analysis suggests a major 
merger between two sub-clumps of similar mass located within the 
central 300 kpc. 

X-ray observations further strengthen our view of a complex 
structure of the inner regions, revealing a hot front in the South- 
Eastern area, most probably a remnant of the occurred merger. This 
merging framework, which arises naturally from our strong lensing 
model, can also reconcile the observed discrepancy between dy- 
namical mass estimates and X-ray, lensing and SZE ones, which 
instead give consistent results. Agreement between X-ray and lens- 
ing mass estimates further indicates that the gas might have had 
time to virialize. Spectroscopic measurements additionally suggest 
that the merger is taking place in the plane of the sky. 

Whereas the presence of a merger is confirmed on several 
grounds, its properties are though still unclear. The detailed fea- 
tures of our model are strictly related to the selection of the mem- 
bers of the multiple image system, in particular to the inclusion of 
the arc AC. Despite all candidates should be confirmed spectro- 
scopically, a photometric analysis suggests that A4 and A5 almost 



undoubtedly belong to the same system. Remarkably, based on the 
shapes and orientations of A4 and A5 alone, without any constraint 
on AC, we can exclude the presence of a single mass component. A 
spectroscopic confirmation of further arc candidates is the required 
step to confirm and accurately define the merging scenario and to 
provide a final, more accurate description of the dynamical state of 
RXJ1347.5-1145. 
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